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Abstract 

The examples of solutions of the Equations of Nonstationary Transonic Gas Flows are 
considered. Their properties are discussed. 



1 Introduction 

Two-dimensional equation of Nonstationary Transonic Gas Flow has the form 

d 2 ( d \ d 2 d 2 

-f(x, y, z) + — f(x, y, z) -£-zf(x, y, z) - ^f(x, y, z) = 0, (1) 



oxoz \ox J ox z oy z 

where variable z is considered as the time-variable . 

Three-dimensional generalization of this equation is defined by the equation 

d 2 id \ d 2 d 2 

2 6^ /(x ' V ' Z ' S) + {dx'^*' V ' Z ' S) ) 6^ /(x ' V ' Z ' S) ~ 6V ^' Z ' S) ~ 

~f(x,y,z,s) = 0. (2) 

where variable z here is considered as the time-variable. 

The solutions of these equations and a corresponding bibliography have been considered re- 
cently in [1] 

In this article we apply the method of solution of the p.d.e.'s described first in [2] and developed 
then in [3], [4]. 

This method allow us to construct particular solutions of the partial nonlinear differential 
equation 

F(x, y, Z, f x , fy, f z , fxx-, fxyi fxz-, fyy-i fyz-, fxxxi fxyyi fxxyi ••) 0. (3) 
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with the help of transformation of the function and variables. 

Essence of method consists in a following presentation of the functions and variables 

f(x,y,z,s) -> u(x,t,z,s), y -> v(x,t,z,s), f x ^u x -—v x , f s -*u s -—v s , 

v t v t 

f - „ - % f - ^ f - 3& f - ^ ~ m 

where variable t is considered as parameter. 
Remark that conditions of the type 

fxy fyxi fxz fzxi fxs fsx--- 

are fulfilled at the such type of presentation. 

In result instead of equation (3) one get the relation between the new variables u(x,t, z) and 
v(x,t,z) and their partial derivatives 

*(u, v, u x , u z , U t , u s , v x , v z , v t , v s ...) = 0. (5) 

This relation coincides with initial p.d.e at the condition v(x, t, z,s) — t and takes more general 
form after presentation of the functions u, v in form u(x, t, z, , s) = F{uj,uj t ...) and v(x,t,z,s) = 
Q{uj,oj t ...) with some function uj{x,t,z,s) . 

Example. 

The equation of Riemann wave 

d d 

faf( x > y) + f( x > y^f( x > y^ = 

after (u, v ^transformation takes the form 

( d \ d ( d \ d d 

The substitution here of the expressions 

d d 
v(x, t) = t-Qj_u(x, t) - u(x, t), u(x, t) = -qW{x, t) 

give us the linear equation 

d . . d . . 

dx^' + di u ( x,t > = 

with general solution 

u(x,t) = _Fl(t-x), 

where „F1 (t — x) is arbitrary function. 

Choice of the function _F1 (t — x) and elimination of the parameter t from the relations 

y - tD(_Fl)(t - x) + _Fl(t - x) = 0, f(x,y)-D(_Fl)(t-x) =0 

lead to the function f(x,y) satisfying the Riemann wave equation. 
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2 Two-dimensional case 

The equation (1) after applying (u,v)- transformation with conditions 

d 

u(x, t, z) = t—u>(x, t, z) — u(x, t, z) 

and 

, . d 
v(x,t,z) = —u{x,t,z) 

takes the form 

( d 2 \ ( d \ d 2 ( d 2 \ d 2 

W 2 ^ M) ^ *' z) M) " 1 " 2 dT 2 ^ *' z) dxTz^ *' z)+ 



( d 2 \ d 2 ( d 2 \ 2 d 

+2 [m; w{x > *' z) ) dwz^ *' z) - 1^^' t? 2) J ^ *' z) = °" (6) 

Its solution of the form 

u(x, t, z) = A(t, z) - xB(t) 
lead to the equation on the function A(t, z) 

' d „ \ d 2 . ( d 1 



with solution 



A(t,z)=_F2(t) + _Fl(z)+ (-1/2 j(jB{t)^ eft + 1/4 (S(t)) 2 j 



where _F2(t), B(t), -Fl(z) are arbitrary functions. 
In result we find that the function 

u(x, t, z) = _F2(t) + -Fl (z) +1-1/ 



2 j(j t B (t)j dt + l/A(B(t))Az-xB(t) 



is the solution of the equation (6). 

After the choice of the functions _F2(t), B{t) and elimination of the parameter t from the 
relations 

( d \ d 

f ~ I * *> z ) ~ w ( x ' *> z ) j = f ~~ ^ t,z) = 

one gets the solution of the equation (1). 
Let us consider some examples. 
In the case 

_F2(t) = 0, B(t) = r\ 

we find the relations 

4 /(x, y, z)t 2 - 4/3 ^t 5 - 8 xt + 3 z + 4 _Fi (^)t 2 = 0, 
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and 

2 yt 3 - zt 5 - 2 xt + z = 0. 

Elimination of the parameter t from these relations in the case _Fi (z) = lead to the solution 
f(x,y, z) of the equation (1) satisfying the algebraic equation 

248832 (f(x, y, z)f z 2 - 248832 x 2 (f(x, y, z)f z + (-1451520 z 2 xy - 221184 y 3 z) (f(x, y, z)f + 

+ (-216000 z A x + 475200 z 3 y 2 + 1327104 yzx 3 + 221184 y 3 x 2 ) (f(x, y, z)f + 
+ (90000 yz 5 + 995328 y 4 zx + 614400 z 3 x 3 + 1290240 y 2 z 2 x 2 ) f(x, y, z)- 

-518400 y 3 z 3 x - 393216 x 5 z 2 - 1179648 x 4 y 2 z + 3125 z 7 - 373248 yV- 
-192000 z 4 x 2 y - 884736 yV = 0. 



3 Three-dimensional generalization 

In a three dimensional case the equation of Nonstationary Transonic Gas Flow takes the form 

d 2 d 2 

f(x,y,z,s) - —f(x,y,z,s) = 0. (7) 



dy 2 ds 2 

Recall that the variable z in this equation play the role of a time-variable. 
After application of u, ^-transformation of the form 

d d 
u(x, t, z, s) = t—tu(x, t, z, s) — iv(x, t, z, s), v(x, t, z, s) = —uj(x, t, z, s) 

we find from (7) the equation 

( d 2 , A d 2 , . n ( d 2 , A d 2 , 

( dj? u ( x > *> z > s ) ) t; z ' s ) - 2 [dj? u ( x > tj z > s) ) *' z ' s)+ 

( d 2 \ ( d \ d 2 

+ \QpW(x,t,Z,s) I \—u(x,t,Z,s) I _ W (a;,i,2!,s)- 

/ <9 2 \ 2 d ( d 2 \ d 2 

~ \m-x^ *' *' s) J ^ ( *' *' *' s) + 2 *' *' s) ) dWz^ *' *' s) " 

From here in the case 

a;(x, i, -2, s) = A(t, s) + k (x + z)t 
one gets the Monge- Ampere equation on the function A(t, s) 

i -(I^H^ (f, ^^ (i ' s) ) 2=o ' (9) 



It is possible to show that the equation (9) can be integrated with the help of corresponding 
(u, v ^transformation. 

Its solutions are dependent from solutions of the linear Laplace equation. 
In fact, the equation 

' if A <)- / d 2 ^ 2 



;f(x,y))—f(x,y)-l—-f(x,y)\ -1 = 



\dx 2 ' J dy 2 ' \dxdy' 
after (u, v ^transformation with 

u(x, t) = tuj t — u>, v(x, t) = UJ t 

takes the form of linear Laplace equation 

d 2 d 2 
—u;(x,t) + —uj(x,t) = 

and its particular solutions after elimination of parameter t give us the solutions of the Monge- 
Ampere equation (9). 

The substitution of another form 

u(x, t, z, s) = A(t, z) + (s + x) t 

into the equation (8) lead to the equation on the function A(t, z) 

having the general solution 

A(t, z) = _F2(t) + _F1 (z) +tz + l/At 2 z 

where _F2(t), -Fl(z) are arbitrary functions. 

The choice of the functions _F2(t), -Fl{z) allow us to construct solutions of initial equation. 
For example in the case 

_F0(*)=r\ _Fl(z) = 
elimination of the parameter t from the relations 

4 f(x, y, z, s)t - t 3 z + 8 = 0, 2yt 2 -2 t 2 z - t 3 z - 2 t 2 s - 2 t 2 x + 2 = 

give us the solution of the equation (7) satisfying the algebraic equation 

16 (f(x,y,z,s)f z+ 

+ (32 zy -32zx- 32 zs - 32 sx + 32 sy - 16 y 2 - 16 z 2 - 16 s 2 - 16 x 2 + 32 yx) (f(x, y, z, s)) 2 + 

+ (72 zy -72zx- 72 z 2 - 72 zs) f(x, y, z, s) + 192 zs 2 + 64 s 3 - 64 y 3 - 
-384 zyx + 384 zsx - 384 zys + 27 z 2 + 192 sx 2 + 64 a; 3 + 64 z 3 + 192 z 2 s + 192 z 2 x + 192 zy 2 + 
+ 192 sy 2 + 192 y 2 x - 192 yz 2 - 192 ys 2 - 192 yx 2 + 192 zx 2 + 192 s 2 x - 384 ysx = 0. 



4 The case of axisymmetric equation 

The equation (7) 

d 2 ( d \ d 2 

2 g^f( x ' y> z > s ) + ( faf( x > y> z > s ) ) d^f^' y > z > s ) ~ A f( x > y > z > ^ = °' 

where 

A — — 

dz 2 ds 2 ' 

in polar coordinates s = rcos(0), z = rsin(0) takes the form 

2 &y^' «■ 2 > + «■ *>) *) " «■ 2 > " = < 10 > 

where we replace variable r on variable z 
After applying (u, v )-transformation with 

d d 
v(x, t, z) = t—u(x, t, z) — iv(x, t, z), u(x, t, z) = —u(x, t, z) 

the equation (10) takes the form 

( d 2 \ ( d \ d 2 ( d \ d 2 

( d 2 \ d 2 ( d 2 \ d 2 

+ ( *> z ) J ^ z ~q^ {j3 ^ *> z ) ~ t3z [ dtif z u ( x i tj z ^ ) ]w u ( x > tj 



d \ d 2 ( d 2 \ 2 d 

+t 2 z [ —u(x, t, z) \ q^u{x, t, z) + zt 2 \—-u{x, t, z)\ —u(x, t, z)- 



cf 2 
dtdx 

This equation has solution of the form 



2zt 2 — — u(x,t,z) = 0. 



u(x, t, z) = A(x, z)t + B(t) 

where function B(t) is arbitrary and the function A(x,z) satisfies the equation 

( d \ d 2 d d 2 

Z \dx~ A{x ' Z) ) dx^ A{x > Z) ~ d~z Mx > Z) " Z dz^ A{x > Z) = °- (U) 

The solutions of the equation (11) can be obtained with the help of (u, ^-transformation and 
a simplest of them has the form 

C2 

A{x, z) = _C2 - _C1 - _C2 ln(- — ) - x. 

z 
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Using the expression 

u(x, t, z) = A(x, z)t + B(t) 

with a given function A(x, z) and arbitrary function B(t) the solution of the equation (10) can be 
constructed. 

As example, in the case B(t) = ln(t) the elimination of the parameter t from the relations 

C2 

f(x, y, z)t - t_C2 + t_Cl + t_C2 ln(- — ) + tx - 1 = 0, 

z 

y-l + \n(t) = 

lead to the solution of the equation (10) 

f{x, y,z) = - (-e~ y+1 _C2 + e- y+1 _Cl + e' y+1 _C2 ln(^) + e~ v+1 x - lj (e^ +1 ) ~* . 
In the case B(t) = te l by analogy way we find the solution 
/(«, y, z) = -« - _« ln( z_) - x + _£» + 1/2 

+ 1/4. ?/ 



(Lom6ertW(l/2 ^/y))' 
The solution of the equation (11) 



A(x,z) = 2/9^ 
z A 



lead to the function 

w(a;,*,z) = 2/9^ + 5(t) 



z 2 



where B(t) is arbitrary function. 
In the case B(t) = ln(t) we find 



/ Xx,y,z) = 2/9^ + (e-y +i y\ 



z 2 

In the case B(t) = t 2 + t~ 2 we get 



f(x,y,z) = 2/9 



^2y + 2Vy T f^Vy T Ti2 + ^2 y + 2Vy T TJ5y^) x 3 



+2/9 



(y + v/^TT^) v / 2y + 2 v ^ + tl^ 

72 ^ 2 + 18 2 2 y 2 + 18 z 2 y^y 2 + 12 
(y + V^TT2) v /27+2^TTI^ 2 ' 
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